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CONFIDENCE  INTERVALS 


FOR  INDEPENDENT  EXPONENTIAL  SERIES  SYSTEMS 

by 

Gerald  J.  Lieberman  and  Sheldon  M.  Ross 


1 .  Introduction  and  Summary 

Suppose  X^Xg,  ...  ,  Xn  are  independent  identically  distributed 
exponential  random  variables  with  parameter  i.e.,  the  density 

function  is  given  by 


( 0  ,  for  x  <  0 

1  -v 

I  \.e  ,  for  x  >  0  . 


Let  Y^Ygj  •«.  >  Ym  also  be  independent  identically  distributed 
exponential  random  variables  but  with  parameter  \2,  and  assume 
that  the  X's  and  Y's  are  independent.  The  problem  is  to  estimate 

R(t)  =  e  1  . 


The  motivation  behind  this  is  that  if  one  has  a  series  system  with 
two  independent  exponential  components  then  R(t)  represents  the 
reliability  of  the  system  a  time  t,  i.e.,  the  probability  that 
the  system  survives  until  time  t. 

In  Section  2  we  present  some  point  estimates.  In  Section  3 
we  present  a  procedure  for  determining  an  exact  (l -a)  level  low-r 
confidence  bound  for  e  .  In  doing  so  we  also  obtain  an 

interesting  characterization  of  the  minimum  of  two  independent  gamma 
random  variables.  In  Section  k  we  compare  our  procedure  with  others 
presented  in  the  literature.  We  show  that  our  procedure  always  gives 
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smaller  confidence  intervals  than  the  procedure  presented  by  Kraemer 
12]  and  that  it  dominates  (in  the  sense  of  yielding  stochastically 
smaller  confidence  intervals)  a  procedure  due  to  Sarkar  [4].  Other 
procedures  are  also  considered.  Section  5  is  concerned  with  the 
same  probltm  when  the  data  is  of  a  different  form,  including  Type  II 
censoring  at  r  out  of  R.  In  the  final  section  Ve  point  out  that 
our  results  easily  extend  to  the  k-sample  problem. 


2.  Point  Estimators 

n  _  m 

Since  X  =  £  X,/n  and  Y  -  £  Y./m  are  the  "usual"  estimators 
i=l  1  i=l 

of  l/\  and  i/\  respectively,  a  "reasonable"  estimator  (the 

1  2  -(vVt 

maximum  likelihood  estimator)  or  e  is  given  by 


Another  possible  point  estimator  is  the  unique  minimum  variance 
unbiased  estimator  (MVUE) ,  which  can  be  obtained  as  follows: 

Letting 

.J1-  lf  xi>t*Ii>t , 

I  0,  otherwise 


it  follows  that. 


*(xx+x2)t 


E(Z)  =  e 

-(Xj^+XgJt 

i.e.,  Z  is  an  unbiased  estimator  of  e  .  Now,  since  the 

n  m  i 

joint  statistic  £  X.,  £  YJ  is  both  sufficient  and  complete, 
ii=l  i=l  / 

it  fellows  from  the  Rao-Blackwell  (also  referred  to  as  the  Lehmann- 
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(z 

n 

2  x , 

m  \ 

2  h) 

l 

i=l  1 

i=l  1/ 

In  m  i  f  n  m  n 


n  m 

)  >  if  £  X,  <  t  or  £  Y  <  t 

i=l  1  "  i=l 

n  n-1  m  m-1 

2  xr‘  2  v* 

=5 —  ^ «  Ev*.  £  Yi  >  t  > 

I)  X  £  Y  i-1  i=l 

i-1  1  .  L  i-1  1  - 


where  the  last  equation  follows  from  the  well  known  fact  that  given 
n 

J)  X. ,  then  X.  is  distributed  as  the  smallest  of  (n-l)  order 
i-1  1  1 

n  i 

statistics  from  a  uniform  distribution  on  0,  J)  X,  .  (A  similar 

i-1  / 

statement  holds,  of  course,  for  the  Y's.) 

Thus,  either  (l)  or  (2)  may  be  used  as  point  estimators  of 

e  .  However,  neither  estimator  is  useful  in  obtaining  a 

confidence  interval  for  the  reliability.  This  difficulty  arises 

because  the  distributions  of  (l)  and  (2)  depend  on  both  \  and  \ 2 

and  not  only  on  +  \g.  In  other  words,  a  nuisance  parameter  would 

have  to  be  eliminated.  In  the  next  section,  we  present  a  method 

-(\  +Ot 

for  obtaining  an  exact  confidence  interval  for  e 

3.  Confidence  Intervals 

We  shall  use  the  notation 


V  ~  Gamma  (7,0) 

to  indicate  that  the  random  variable  V  is  distributed  as  a  gamma 
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random  variable  with  parameters  y  and  f3,  i.e.,  its  density  function 


is  given  by 


v  >  0 

v  <  0  , 


where  >  0.  Similarly,  the  notation 

v  ~  x2 

r 


will  indicate  that  the  random  variable  V  is  distributed  as  a  chi 
square  random  variable  with  r  degrees  of  freedom. 

Now,  let 


U  =  min 


n 

-\ 

Li 

i-1 


m 


E  Y, 


i=l 


and  let 


K  =  ^largest  j  <  n:  Xi  5  u  }•{  largest  j  <  m:  <  uj 


Note  that 


rain(m,n)  <K<m+n"l* 


We  shall  show  that  given  K  =  k,  then 

U  ~  Gamma  (k,  \^  +  \ 2)  , 
or  equivalently,  that  given  K  =  k,  then 
(3:  2(X1  ~  X^k  . 

Equation  (})  enables  us  to  determine  an  upper  confidence  bound  for 

k 


X^  +  Xg,  with  confidence  coefficient  1-a,  namely 


N 


where  '<  n,  is  such  that 
Ct»2k 

p(4  >  x;U>  - a  • 

From  (4).  we  can  easily  obtain  a  lower  confidence  bound  for  R(t) 

“ ( X^+Xg  ^ 

e  ,  with  confidence  coefficient  1-a,  namely 


(5) 


exp{-Xa,2kt/2u)  * 


The  remainder  of  this  section  will  be  devoted  to  showing  that 
given  K  =  k,U  ~  Gamma  (k,  X^  +  Xg).  To  facilitate  this  we  shall 
slightly  change  our  notation  and  also  make  use  of  the  equivalence 
between  the  partial  sums  of  independent  identically  distributed  expo¬ 
nential  random  variables  and  the  Poisson  process. 

Let  {N^(t),  t  >  0}  and  {Ng(t),  t  >  0)  be  independent  Poisson 
processes  having  respective  rates  X^  and  Xg.  Also,  let  N(t)  = 
N^(t)  +  Ng(t).  It  follows  from  known  results  that  (N(t),  t  >  0} 

is  also  a  Poisson  process  but  with  rate  Xj^  +  Xg.  Now  let  be 

s*t  th 

the  time  between  the  (i-l)  and  i  event  in  the  N(t)  process, 
and  define 

>  th 

1  if  the  i  event  in  the  N(t)  process  canes  from 

the  N, (t)  process, 

Ji  =(  th 

2  if  the  x  event  in  the  N(t)  process  comes  from 

the  Np(t)  process. 

Finally,  consider  a  new  Poisson  process  (N*(t),  t  >  0} ,  inde- 

xy. 

pendent  of  the  above,  and  having  rate  X^  +  Xg.  Let  T*  be  the  i 
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interarri/al  time  for  this  process  and  let  J*,  i  =  1,2,  ...  be  an 
independent  identically  distributed  sequence  of  random  variables, 
independent  of  the  above,  and  such  that 


/  1 
1  with  probability 


1* 


with  probability 


\l  +  ^2  * 


h  +  X2' 


Theorem  1:  The  two  stochastic  processes  i  >  1}  and 

{(T*.»J*)>  i  >  1J  have  the  same  probability  law. 

Proof: 

The  theorem  will  be  proven  by  using  induction  to  show  that 

PITjiV  J1  •  31’  T2  —  a2’  J2  -  V"”’n<-V  Jn“ 

-  P(I*  <  ar  J«  -  Jj,  TJ  <  a2,  J|  *  V-1!  -  V  JS  *  ■>n)  ‘ 


Let 


Now , 


■{ 


2  if  i  =  1  , 
1  if  j  =  2  . 


r°° 

(6)  PIT,  <  a,  J  =  j}  -  /  P{T..  <  a,  J  =  j|Y  =  yj\.  e 

x  ”  1  JO  Jc 


"V 


dy 


where  Yjl  is  defined  to  be  the  first  interarrival  time  of  the  process 

(N.  (t),  t.  >  Oj .  Hence, 

J  „ 


P{TX  <  a,  J1 


r*  -V  V 
=  J }  =  /  (l-e  J  )\,  e  dy 


-00  -\A 

t  ^  A  J  ) 


+  /  (l-e  J  )\.  e  dy 


=  > 

=  P(Tj£  <  a,  J#  =  j}  . 

Note  that  the  last  equality  follows  since  T|  is  just  the  time  of 
the  first  event  of  a  Poisson  process  with  rate  +  \g  and  is 

j  with  probability  - — *-r—  (independently  of  anything  else).  Now 

by  (7)  we  have  that 


P{T!  -  V  Ji  *  Jl'  *  ‘  *  'Tn  -  V  Jn  “  Jn} 


f*l 

’  jl  PfT2  <  a2;  J2  *  J2'-'Tn  1  V  Jn  "  jJTl  “  a'  J1 


J1)X 


-(\1+\2)« 


However,  by  the  lack  of  memory  of  the  Foisson  Process  it  follows 
that  at  time  (which  is  the  first  time  that  an  event  from  either 

N-^  or  Ng  occurs)  both  the  processes  (N^(t),  t  >  0)  and 
(Ng(t),  t  >  0)..— start  over  again,  and  hence  (8)  equals 

r-(\,+\2)a 
P(T1  -  *2'  J1  =  V  •••  '  Tn-1  -  V  Jn-1  =  da 


which  by  our  induction  hypothesis  equals 
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-(\,+\p)a 

da 


Q,  •  E  •  D  • 


Now  consider  the  first  event  in  the  N(t)  process  such  that  to  the 
left  of  (and  including)  that  event  there  are  either  n  events  having 
their  J's  equal  1  or  m  events  having  their  J's  equal  2.  Suppose 
this  is  the  Kth  event  and  let  U  *  +  •  •  •  +  TR  denote  the  time  of 

this  -vent  (Note  that  these  are  the  same  K  and  U  as  were  defined 
at  the  beginning  of  this  section.) 


FIT*  <  a  , 


1  J2',,,/ln-l  - 


J«V 


Theorem  2:  Given  K  =  k,  U  has  a  gamma  distribution  with  parameters 

k  and  .  Tnat  is,  its  conditional  density  is  given  by 

1 


.ft.  i  , 

-  X  A  k-1  -\U 

i K=k*'u  =  (k-l)r  u  e 


Pr  oof 

p.  fine  K*  ana  U*  in  an  equivalent  manner  for  the  N*(t) 
p,  ot'csu  Since  the  values  are  determined  completely  at-  random 

for  the  N>  process  the  result  follows  for  U*  and  K*.  The  result 
must  t. neri  follow  for  U  and  K  by  Theorem  1. 

Q.E.D. 
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Corollary  3;  Let  X  ~  Gamma  (n,\^)  and  Y  ~  Gamma  (m,\g),  where 
the  X's  and  Y's  are  independent.  Then 


where 


P{K  *  k}  =  | 
1 


min  (X,Y)  ~  Gamma  (K>\^  +  \g) 


k-l|  I  *1  in|  ^2  lk"n  +  fh-lj  \  k"m  __ 
n-lll^+^J  (Xj+XgJ  lm-11  \1+\2J  \x' 


\  |k-m/  \2  im 


for  min  (m,n)  <k<m+n-l. 

The  correct  way  of  interpreting  Corollary  3  is  that  conditional  on 
K  =  k,  min(X/Y)  ~  Gamma  (k,^  +  kgb  i-e.,  letting  rk^(a)  be 
the  probability  that  a  gamma  random  variable  with  parameters  (k,\) 


is  greater  than  a,  then 


m+n-1 


P(min  (x,Y)  >  a}  ,  £  rk  v  ( 

k=min(m,n)  *'V*2 


a)p(K  »  k} 


Proof: 


The  proof  is  immediate. 


Comparison  with  Other  Procedures 


Sarkar  [4]  considered  the  case  ra  =  n  and  obtained  an  exact 

-f  l+^Jt 

confidence  bound  for  \1  +  \g  (and  hence  for  e  )  by  the 

following  method.  He  defined  =  min(X^Y^),  i  =  1,  ...  ,  n,  and 

based  his  confidence  bound  on  the  fact  that  the  are  independent 

identically  distributed  exponential  random  variables  with  mean 

n 

l/(\l+\p).  Hence.,  £  Z i  ~  Gamma  (n^+Xg).  His  upper  confidence 
bound  for  \1  +  \2>  with  confidence  coefficient  1-a,  is  thus 
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i,2n'2  |  Zi  • 

The  weakness  of  Sarkar's  procedure  is  that  he  is  throwing  away  a  good 
deal  of  information.  For  instance,  if  <  Y^  then  he  only  uses 
and  makes  no  use  of  Y^.  However,  given  that  X1  <  Y1  it  follows 
(by  the  lack  of  memory  of  the  exponential)  that  the  excess,  Y^  -  X^ 
is  exponentially  distributed  with  mean  l/\g.  Hence  by  using  the 
"excesses"  one  can  generate  more  X's  and  Y's.  This  is  equivalent 
to  what  the  procedure  of  Section  3  does.  Let  us  illustrate  by  an 
example . 

Example :  Independent  random  samples  of  size  4  are  obtained  on  each 
of  two  components,  X  and  Y.  Assume  that  the  X's  are  exponentially 
distributed  with  parameter  \1  and  the  Y's  are  exponentially  dis¬ 
tributed  with  parameter  kg.  Suppose  the  data  are  as  follows: 

(x1,x2,x3,x4)  =  (2, 7.1, 6, 3) 

(WVV  =  (5A,2,3.1)  . 

Sarkar  would  thus  generate  a  random  sample  of  size  4,  namely,  (2,4, 2, 3) 
from  an  exponential  distribution  having  mean  l/^+Xg).  Our  procedure 
(described  in  Section  3)  is  equivalent  to  the  following:  Determine 
the  minimum  of  X1  and  Y^,  thereby  obtaining  the  value  of  2. 
however,  note  that,  the  excess,  Y^  -  X^  =  3>  yields  another  Y  data 
point.  Hence,  after  generating  the  value  2  we  are  left  with  the 
"random  samples" 

X:  7-1,  6,  3 
Y:  3,  4,  2,  3.1  • 
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Again,  we  determine  the  minimum  of  the  first  X  and  the  first  Y, 
i.e.,  min(7.1,3)  =  3,  thereby  obtaining  the  "outcome"  3  and  note 
that  the  excess,  7*1  -  3  =  4.1,  yields  another  X  data  point. 
Having  generated  the  value  J,  we  are  left  with  the  "random  samples" 

X:  4.1,  6,  3 
Y:  4,  2,  3-1  • 

We  determine  the  minimum  of  the  first  X  and  the  first  Y,  i.e., 
min  (4.1,4)  =  4,  thereby  obtaining  the  "outcome"  4  and  note  that  the 
excess,  4.1  -  4  =  0.1,  yields  another  X  point.  Having  generated 
the  value  4,  we  are  left  with  the  "random  samples" 

X:  0.1,  6,  3 
Yi  2,  3.1 


The  minimum  of  the  first  X  and  the  first  Y  is  0.1,  thereby  obtaining 
the  "outcome"  0.1.  The  excess,  2  -  0.1  *  1.9#  yields  another  Y 
point.  Having  generated  the  value  0.1,  we  are  left  with  the  "random 
samples " 

X:  6,  3 
Y:  1.9#  5-1  • 

Since  min(6,1.9)  *  J..  9#  we  obtain  the  "outcome"  1.9,  and  the  excess 
of  4.1  yields  another  X  point.  Having  generated  the  value  1.9,  we 
are  left,  with  the  "random  samples" 


X:  4.1,  3 


Y:  3.1 
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Since  min(U.  1,3.1)  =.  3.1,  we  obtain  the  "outcome"  3-1;  and  the  excess 
of  1  yields  another  X  point.  Having  generated  the  value  3.1,  we 
are  left  with  the  "random  samples" 


X:  1,  3 
Y: 

Since  there  are  no  more  Y  data  points  we  must  stop  the  procedure. 

If  we  treat  the  "outcomes"  as  the  outcomes  of  a  random  sample  of 
size  6  -  namely  2,  3*  **,  0.1,  1.9*  3-1  -  from  an  exponential  distri¬ 
bution  with  mean  l/fX^+Xg),  then  the  sum  should  behave  as  a  gamma 

distributed  random  variable  with  parameters  6  and  X^  +  Xg.  However, 

In  m  . 

X^,  Yj 
i*l  i*l  / 

defined  in  Section  3,  and  6  is  just  the  value  taken  on  by  the  random 
variable  K  a^so  defin-d  in  Section  3»  From  the  results  of  Section  3 
these  generated  values  can  then,  in  fact,  be  considered  as  the  outcomes 
of  a  random  sample  from  an  exponential  distribution  with  mean 
i/(X,+\~).  As  indicated  earlier  our  confidence  bound  would  be  arrived 
at  by  using  the  fact  tnat  U  ~  lamir.a  (K  =  6,X1+Xg)*  Hence,  our  pro- 
<  -..dure  yields  stochastically  smaller  confidence  intervals  than  does 
Jarkar's  since  w.  ar-;  generating  a  larger  random  sample  (i.e., 
n  <  k  <  : 

The  above  example  also  points  out  some  of  the  weaknesses  of  our 
procedure.  First,  since  we  had  no  Y  points  to  compare  them  with, 
we  had  to  throw  out  the  final  two  X  points.  Also  K,  the  number 
of  dat.a  points  used,  is  a  function  of  the  ordering  of  the  X's  and 
Y's.  Thus,  for  instance,  random  samples  (X^,Xg,X^) ,  (Y^,Yg,Y^) 
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might  product  a  different  value  of  K  than  would  random  samples 
(x2,xrx3),  (Y2,yi,Y5). 

Another  procedure  is  the  one  due  to  Kraemer  [2],  Kraemer  uses 
n 

the  fact  that  ^  X.  ~  Gamma  (n,\.)  (or  equivalently  that 

1  1  1 
n  „  m 

2\1  ~  Xgn)  and  Yi  Gamma  (m,\2).  Thus,  since  the  X's 

and  Y's  are  independent  it  follows  that 


‘  n  m 

H  £ 

.  1  1 


~  X! 


'2(n+m)  1 


have  that 


in  m  » 

Again,  letting  U  «  min|2  Xi>  £  YJ  we 

In  m  I 

xi +  x2  r  Yij  • 


However 


[’  n  m  I 

hj 


~  X2(n+m)  S0  that  an  upper  confi(ience 


bound  for  with  confidence  coefficient  at  least  1-a,  is 

given  by 


^  Xa;2(n+ra),/2U  * 

Equation  (9)  expresses  the  confidence  bound  for  Kraemer's  procedure. 
However,  our  procedure  (equation  (4))  yields  a  smaller  confidence 
interval  since  K<  2(n+m).  This  is  intuitively  clear  since  Kramer 
bases  her  estimate  on  U  and  gets  a  bound  on  its  distribution  while 
we  base  our  estimate  on  both  U  and  K  and  use  the  exact  joint 
distribution  of  U  and  K. 

Another  procedure  has  been  suggested  by  Lentner  and  Buehler  t .3 J . 
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•» 

1 


n  m 

Their  procedure  is  based  on  the  fact  that  given  2/  ^  ~  /)  >  then 

i  i 

n 

the  distribution  of  £  X.  depends  on  X  and  only  through 

1  1  id 

Xi  +  \p.  They  then  use:  the  Lehmann-Scheffe  theory  to  obtain  a  con- 
n  m 

ditional  (on  2/  X.  ~  Ti  uniformly  most  powerful  unbiased  test 

1  1  .1  1 

for  the  hypothesis  V  *•  \g  s.  c.  By  inverting  this  test  they  obtain 
an  upp^r  confidence  bound  for  with  confidence  coefficient 

1-a. 

n  n  m 

Since  the  conditional  distribution  £  X.  given  £  X.  -  £  Y. 

1  1  1  1  1 

is  somewhat  complicated  it  is  not  theoretically  clear  how  good  the 

Lentner-Beuhler  procedure  is-  It  should  also  be  pointed  out  that 

their  methods  are  computationally  difficult- 

other  relevant,  literature  that  should  be  mentioned  is  a  procedure 

given  by  Grubbs  [1]  which  obtains  fiducial  intervals  for  the  reliability. 

Anothei  possible  approach  to  the  problem  would  be  to  assume  prior 

aistributionson  and  \0.  It  would  then  be  straightforward  to 

X  d 

obtain  the  posterior  distrioution  of  ,\^  +  \ ^  and  Bayesian  probability 
intervals  can  then  t  obtained. 


J. 

I 

l 

1 

I 

i 


i 

1 

1 
1 

2 

T 

** 


5  Otter  Sampling  Srhvmcs 

common  ryp-  of  sampling  used  in  practice  is  type  II  censoring 
at  r  out  of  ft  (r-  <  ft)  Type  II  censoring  at  r  out  of  R  refers 

to  the  situation  where  ft  it  ms  are  tested  simultaneously  and  the 
test  is  terminated  wren  tne  first  r  out  of  R  have  failed.  The 
time  of  failure  of  each  item  is  recorded,  This  corresponds  to 
recording  tr>  order  statistics,  i.e„,  X^^X(gi,  ...  ,  X^j .  In 
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the  notation  of  this  paper  put  N  units  of  the  first  component  (x) 

XL 

on  test  and  terminate  the  test  when  the  n  failure  occurs.  Similarly 
put  M  units  of  the  second  component  (y)  on  test  and  terminate  the 
test  when  the  m  11  failure  occurs.  Let  and  Y^j  be  the 

respective  order  statistics  (i  =  1,2,  ...  ,  n  and  j  =  1,2,  ...  ,  m). 
Denote  by 


6i  *  (n-i+l)(X(i)  -  X(i-1)),  i  -  1,  2,  ...  ,  n  (X(Q)  s  0) 
and 

j  j  <)  ■  2,  •••  >  *  (Y(0)  s  0)  • 


It  is  well  known  that  the  6's  and  t's  are  independent  and  identi¬ 
cally  distributed  random  variables  with  respective  parameters 
and  \g.  Hence,  the  procedure  of  Section  3  is  applicable. 

Another  type  of  sampling  can  be  described  as  follows.  For  the 
first  component  (x)  put  one  item  on  test.  When  it  fails  replace  it. 
When  this  second  item  fails  replace  it.  A  similar  procedure  is 
followed  for  the  second  component  (Y).  Terminate  the  test  when 
there  have  been  either  n  failures  of  component  1  (X's)  or  m 

failures  of  component  2  (Y's).  Follow  the  techniques  outlined  in 

-(\  +\g)t 

Section  3  to  determine  a  confidence  bound  for  e  .  Note 

that  in  this  case  our  procedure  uses  all  the  available  information. 

The  data  can  be  represented  by  X, ,  . .  •  ,  X  and  Y. ,  ...  ,  Y 

x  r  x  s 

where  either  r  =  n  or  s  *  m.  Consider  the  case  r  =  n,  then  the 
likelihood  is  given  by 


L(x1,...,xn,y1,...,yg)  =  e 


’*2 
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where  the  first  term  arises  since  the  next  y  value  -  namely  yfl+^  * 

n  s 

must  be  greater  than  2/  x<  "  £  >V  Hence, 

11 

n 

L(xl'  *  Vyl*  *■'  '  ys^  “  XlX2e 


Letting  \  *  \  +  \g  and  6  >  we  have 

n 

'Ij  x 

n  s  1  A 

L(x^,  •  •  •  ,  xn,yl*  *  ys^  *  ®  (^“®)  * 
or 

log  l(x1#  . . .  ,  xn,y1#  •  •  •  ,  yg)  ==  n  log  0  +  s  log  (\-0)  -  \ 


Taking  the  partial  derivatives  and  setting  them  equal  to  zero  leads 
to 

K 

4: log  L  =  rH  -  £  xi " 0 

and 

d  ,  .  n  s 

3e  log  L  =  e  ’  r^e  =  0  * 

Solving  these  equations  simultaneously  results  in  the  maximum  likeli¬ 
hood  estimate  \  of  \,  i.e., 


A  s  +  n 
\  =  - 


a 

1 


Also,  in  the  case  where  s  =  m  we  have  that 


r  +  m 
m 
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Suppose  X,  ,Xo, . .  •  ,x‘  are  independent  identically  distributed  exponential 
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random  variables  with  parameter  \^.  Let  Y^,Yg, . . .  ,Y;k  also  be  independent  iden¬ 
tically  distributed  exponential  random  variables  with  parameter  \g,  and  assume 
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exponr’vtial  components  then  R(t)  represents  the  reliability  of  the  system  at 
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A  procedure  for  determining  an  exact  (l<x)  level  lower  confidence  bound 
for  R(t)  is  presented.  In  doing  so  an  interesting  characterization  of  the 
minimum  of  two  independent  ganma  random  variables  is  obtained.  The  suggested 
procedure  is  then  compared  with  others  presented  in  the  literature. 
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